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Abstract 

We develop a computational model of quantum turbulence decay employing a kinematic pre- 
scription for the normal fluid. We find that after an initial transient, the length of the vortex tangle 
L decreases and for large times obeys the scaling law L ~ t~ 0A5 . The average magnitude (along 
the quantized vortices) of the superfluid and line-vortex velocity are close and differ significantly 
from the average magnitude of the normal fluid velocity. 

PACS numbers: 67.40.Vs, 47.27.Ak, 47.27.Gs 
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In quantum turbulence physics a tangle of quantized vortices interacts via mutual 
friction forces with thermal excitations (normal fluid) of the superfluid ground state. An 
elementary research program for this multifaceted problem investigates idealized flows 
characterized by symmetries like homogeneity in space and/or time, as well as, isotropy. 
At first, phenomenological issues like scalings of energy spectra , energy decay |J| and 
structure functions should be established. This Letter contributes to this research program 
by employing a mathematical model of decaying quantum turbulence under conditions 



resembling the experiment of 
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and solving it with numerical and computational methods. 



Our model consists of a dynamic equation describing the superfluid vortices and a kine- 
matic prescription for the turbulent velocity field. In particular, if S(£,t) is the three 



dimensional representation of the vortex tangle then its motion obeys the equation 

— = V t = hV s + h x S'x(V n -V a )- 



.1: 



h xx S'x(S'xV n ) (1) 



where the superfluid velocity V s is given by the Biot-Savart integral: 

k f (S — x) x dS 



V '( x ) = V , ' l3 , (2) 
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where t is time, x is space, k is the quantum of circulation, V n is the velocity of the normal 
fluid, £ is the arclength along the loops, S' = p^rr is the unit tangent vector while h, h x and 
h xx are constants related to mutual friction physics. 

At every instant, the normal velocity is decomposed into a mean value and a fluctuation 
V n = (V n ) + u n with the fluctuation u n defined by the following function Q]: 

M 

u n = ^2\A m x k m cos(k m ■ x + uj m t) + 

m=l 

B m x k m sin(k m ■ x + u m t)], (3) 

where M is the number of wavemodes the sum of which constitutes the velocity field. A m , 
B m are vectors with random orientation and magnitude | A m | 2 = |-B m | 2 = (2/3)E n (k m )Ak m 
with E n (k m ) the normal fluid energy spectrum at wavenumber k m . In addition, k m is a unit 
vector normal to both A m and B m and k m = k m k m . The frequencies uj m = yl^EdJz^) 
correspond to the physical notion of "eddy turnover time" . As required, u n is incompressible 
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FIG. 1: Vortex tangle at initial, maximum length and stoppage times. To serve clarity, only one 
eighth of the computational box is shown. 

by construction. The energy spectrum E n (k m ) has two parts: A high wavenumber part 
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("inertial range"): E n (k m ) = CK^ 3 k m 3 , where Ck = 1.5 is the Kolmogorov constant and e 
is the rate of energy dissipation and a low wavenumber part ("large eddies"): E n {k m ) = Ak^, 
where A is defined by matching the two spectra at k e the "integral length scale" wavenumber. 
According to the permanence of large eddies hypothesis, A remains constant during the 
decay of turbulence. We observe that in the present definition of u n there is no coupling 
between different modes and so there can be no energy flux in wavenumber space. This is in 
contrast to what happens in Navier-Stokes turbulence. In addition, the model is insufficient 
for capturing fine normal fluid effects like intermittency. On the other hand, it reproduces 
adequately a number of turbulence phenomenologies that relate to even order statistics like 
for example the Kolmogorov scaling or the Lagrangian flatness factor Moreover, since 
the model was devised having homogeneous, isotropic turbulence in mind, it is suitable for 
the description of the grid normal fluid turbulence in the jj| experiment. 
We use the symbol l e for the integral scale of turbulence (peak of the spectrum) and the 
symbol u' n for the turbulence intensity. It is 3u' n 2 = Y^=i( u n u n) = where E is the 
kinetic energy of turbulent fluctuations. Knowledge of l e and u' n at each time step allows 
the construction of the normal velocity field. In particular, we can calculate the turbulent 
Reynolds Ret = u' n l e /u (where v is the coefficient of viscosity), the energy dissipation rate 



e = u' n 6 jl e and the Kolmogorov scale r\ = l e /Re^ 4 . We employ the decay model of |4( in 
order to calculate l e and u' n as functions of time. According to this model, there are two 
periods of turbulence decay: During the first, l e is smaller than Z& (the box size). During 
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the second (which starts at t s = — ^s-y ~a(^ 2 — l% 2 ) with l e0 being the integral scale at 

5(27r)2 V 

t — 0), l e remains constant and equal to The change in l e before its saturation is given by: 

/ \ 2/5 

l e (t) = 2tt( (5/ll)y^4(t + t )J , with t = -^-^\l^fl% 2 . For times smaller than the l e 

saturation time the kinetic energy of turbulence E is calculated from: E(t) = Eq(1 + ^) , 
where Eq = (9/6)(27r) 3 A/Z 3 . We have found that the prefactor 9/6 is necessary in order not 
to have a discontinuity in normal fluid energy at t s . This condition is not satisfied by the 



prefactor 11/6 in formula (4) of |4|. For post-saturation times it is: E(t) = 2 , 2 ^ (t + t + 
ti) -2 , with t\ = (4/5)(27r) _5//2 C^ 2 y4 _1//2 /^ 2 . The constants to and t\ define the virtual origin 
time t vo — — (t — ti). As defined in jgj t vo is the time when (supposedly) the turbulence 
has infinite energy concentrated on an integral length scale of infinite wavenumber. In 
this interpretation as the energy decays the energy containing wavenumber moves towards 
smaller values. Our initial conditions correspond to an intermediate turbulence state in 
the decay process. We employ periodic boundary conditions for the superfluid tangle by 
introducing image vortices. The normal flow is periodic by construction. 
The working fluid is 4 He — II and so the quantum of circulation has the value k = 9.97 • 
10~ 4 cm 2 / s. The calculation is done at T = 1.3K (compared with T = 1.5K in 0) for which 
the other parameters of the problem have the values: v = 23.30 • 10 _4 cm 2 /s, h = 0.978, 
h x = 4.0937 • 10 -2 and h xx = 2.175 • 10~ 2 . In addition, we have Re t = 5 • 10 3 and the initial 
peak of the spectrum is located at Z e o = 0.0161cm which corresponds to k e0 w 62cm" 1 . For 
comparison, = 0.1cm and = 10cm -1 . We mesh the line vortices with discretization 
length Aa; = 4/84 = 1.19 • 10 _3 cm. We use the same distance to define the smallest 
resolvable wavelength in the normal fluid turbulence model (eq. (3)), l co = 2.381 • 10~ 3 cm 
which corresponds to wavenumber k co = 420cm -1 . Using the equations of the model we can 
calculate the saturation time t s = 0.2048 ■ 10 -2 s, as well as, the time t vr at which 7] will 
become equal to the smallest resolvable scale l co : t vr ~ 0.025s. Beyond this time the line 
vortices could develop structure at space scales smaller than the smallest normal turbulence 
wavelength. At stoppage time t e = 0.389s, the Kolmogorov scale is equal to 5.26-10~ 3 cm and 
therefore it is greater than the smallest resolvable wavelength equal to 2- Ax = 2.38-10 _3 cm. 
The stoppage time Kolmogorov scale corresponds to wavenumber k ve = 190cm _1 . The time 
step is chosen in order to ensure that none of the (resolvable by the numerical grid) Kelvin 
waves propagates more than one discretization vortex segment within one calculation step. 
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The typical time step for this is At ~ 0.3 • 10 _3 s. Initially, the tangle consists of 114 vortex 
rings of random orientation and its total length is Lq = 14.06cm. The choice of random 
initial conditions is justified for two reasons: (a) there is no experimental information about 
the actual initial tangle configuration which could be employed, (b) since the experimental 
scalings are reproducable without an explicit control over the geometry of the initial vortices 
it must be the case that the scaling phenomenology does not depend on the latter geometry. 
The second point agrees with our conception of turbulence as a statistical flow state that 
can be achieved from a variety of initial conditions. This is equivalent to our understanding 
of turbulence properties as idiosyncratic of the differential equations governing the system 
and not of the initial conditions. Different initial conditions are driven by the system to 
generic (reproducable) turbulence scalings. 

Other useful quantities one can calculate are the average values of the velocity magnitudes 
| I'll, \V n \, \V a \ and |cos(0)| = \S' ■ (V n — Vi)\/\ V n — Vj| along the vortex filaments at various 
times during the system's evolution. These averages are taken by sampling the quantities 
of interest at each discretization node, and subsequently forming their arithmetic mean. 
Angle 9 is an important quantity in the physics of the mutual friction force per unit length: 
/ = p s nd xx S' x (S' x (V n - Vi)) - p s nd x S' x (V n - V t ). Here d x = -2.045 ■ 10~ 2 and 
d xx = 4.270 • 10~ 2 are nondimensional coefficients and p s = 138.6 • 10~ 3 g/cm~ 3 is the 
superfluid density. 



The results {(V n ) = 0) show that in accordance to an instability discovered by Cheng 
et al and elaborated mathematically by Glaberson et al ^^|, the normal flow excites 
Kelvin waves on the filaments (Fig^). The length of the tangle reaches a maximum of 
Lmax = 47.83cm at t — 0.0019s and subsequently decreases. At maximum length the 
vortex line density is A = L/V ~ 0.5 • 10 5 compared with A pa 2 • 10 5 for the smallest A 
run in J|. In stating this, we made use of the relation uj(t) = kA in order to deduce A 
from their superfluid vorticity (co(t)) data. Refering again to the smallest A run in 4J, we 
note that A varies there over 2 orders of magnitude while here only by a factor of 5. As 
required (Fig|21 right), the decay of turbulent normal fluid energy obeys the two previously 
mentioned temporal scaling laws. Also demonstrated in the same figure (left) are the two 
spatial scaling regimes in the E n {k) spectrum before t s , as well as, the disappearance of 
the large eddies scaling regime for subsequent times. One can ask a two fold question: (a) 
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FIG. 2: Left: normal fluid spectra E n (k) at t = 0, t = t s , t = t^ r and t = t e . Right: the evolution 
of tangle length L and normal fluid turbulent energy E n . 

why the vortex length decreases after a rapid transient and (b) why its observed temporal 

scaling at large times, L « t~ 0A5 , differs from the L « t~ 1,5 one of y? 

Possible reasons for the latter might be the inadequacy of the employed turbulence model or 

the shorter decay time span of 2 decades in the calculation compared to 3 in the experiment. 

In this milieu, an important question has to do with the meaning of the reported tangle 

I I 

lengths in both theory and experiment in the light of the findings of [llj that the superffuid 
tangle is a fractal. According to [l^ (pg 25), the latter means that, as long as, the yardstick 
for length measurements belongs to the scale range within which the tangle satisfies a fractal 
scaling, different yardstick lengths will give a different length for the tangle. When it comes 
to calculation, this implies that better resolved fractal tangles would be measured to posses 
significantly greater lengths when the (different for each resolution) discretization length 
is used as yardstick length. Moreover, one could ask at first whether the second sound 
measurement technique introduces (in fractal turbulent tangles) such a yardstick length 
depending on the second sound wavelength/frequency. In other words, whether there exists 
a certain Kelvin wave frequency above which the experimentally employed second sound 
does not see the variations of the line density in a fractal vortex system. Subsequently, 
whether this possible uncertainty in the length measurement affects the scalings observed 
during the decay of turbulence. These matters are not clear-cut issues and deserve further 
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investigation by expanding (for example) the work of which calculated the motion of 
one roton shot toward a single straight line vortex in the realm of interactions between a 
roton and fractal vortex lines. 

It is useful here to notice that in contrast to classical turbulence a uniform normal velocity 
field is not compatible with the assumption of isotropic superfluid turbulence and it is 
not dynamically irrelevant. These are clearly seen in counterflow quantum turbulence 
calculations [1J| where the superfluid turbulence is due exclusively to a constant (externally 
imposed) velocity field. In these calculations, there is anisotropy in the direction of the 
imposed flow. The above are reminiscent of the role of uniform, imposed, magnetic fields 
in inducing anisotropies and affecting the velocity spectra in hydromagnetic turbulence 
(page 132 in and page 100 in 16\). In order to clarify better this point we have done 
two more calculations, one with stationary turbulence and another with turbulence decay 
but with the addition of a constant normal velocity field in the y direction. We have fixed 
the velocity magnitude = V ny — 7cm/ s so that it is comparable to the grid towing 

velocity (between 5 and 200 cm/ s) in |4j. Although in this way u' n /V ny ~ 100 (at t = 0) one 
observes in FigJHl (right) that with the imposed velocity field the vortex length increases 
at times for which (in the case of purely decaying turbulence) it decreases. Therefore, the 
temporal decay law for the line-vortex length could be affected by small bulk normal fluid 
velocities. We implicitly assume here that the mean velocity profile is stable; this is also 
the case in jl^ . 

To analyze the first leg of the previously posed question, we note that according to the 
analysis of the instability (see also the discussion in j^|) the normal fluid velocity 

fluctuations transfer energy to a Kelvin wave of a particular wavelength (and therefore 
increase its amplitude) only when their component along the direction of motion of the 
wave is both greater in magnitude than the group velocity of the wave and parallel (of 
the same sign) to it. Any normal fluid velocity antiparallel to a vortex wave reduces its 
amplitude. Therefore since initially the turbulence intensity is approximately 100 times 
the group velocity of the fastest (resolvable) Kelvin wave and there are no Kelvin waves 
present, there is unhindered Kelvin wave excitation. This explains the initial rapid increase 
of vortex length. Subsequently, in conjunction with the decay of turbulence intensity, at 
places with wave group velocity (a) larger than the local normal fluid velocity magnitude or 
(b) antiparallel to the normal fluid velocity direction, the wave amplitudes are damped. The 
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FIG. 3: Left: average values of \Vi\, \V n \, \V 8 \ and |cos(0)| along the quantized vortex filaments. 
Right: evolution of tangle length L in three calculations with different normal fluid turbulence 
characteristics. 

turbulence decay factor is necessary since as FigEl (right) shows, in stationary turbulence 
the length keeps increasing rapidly. 

The results for the average values of \Vi\, \V n \, \V S \ and |cos(0)| are presented in Fig 01 
(left). They show that the initial length transient corresponds to rapid (\Vi\) and (|V S |) 
transients. The latter increases towards a maximum coinciding in time with the length 
maximum while the former decreases approaching (|V S |). Notice that when the length 
starts decreasing there is order of magnitude difference between (|V^|) and (|V^|). This is 
particularly true for the average of \V B \ taken over the whole volume of the fluid since this 
was found (using a 84 3 grid) to be an order of magnitude smaller than the average over the 
line vortices. Another observation is that for smaller than (|V S |) (which happens at 

very small normal turbulence energies) S' tends to become normal to V n — Vj,. The results 
bring forward quantum turbulence physics that differ significantly from those proposed 
in where the assumption was made that the superfluid and normal fluid velocities are 
identical. However, one must also bear in mind that the employed turbulence model does 
not have the quality of fully dynamical Navier-Stokes calculations. The latter kind of 
computations would eventually be required in order to verify the present results. 



Overall, it is not possible to argue for the quality of the established mathematical models 
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in quantum turbulence theory, as long as, the available experimental data fail to address 
directly the main variables (fluid velocities, vortex tangle geometry) in these models. This 
is a major obstacle for progress in quantum turbulence research. The present work suggests 
that progress depends on the development of new, more potent experimental methods and 
their combination with fully dynamic mathematical calculations. The latter could identify 
generic and essential phenomenological trends that could be encoded in statistical mechanical 
equations. 
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